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The elliptical stadium is a plane region bounded by a curve constructed by joining two half-ellipses, 
£"""•» i with half axes a > 1 and b — 1, by two parallel segments of equal length 2h. 

V. Donnay [0 proved that if 1 < a < y/2 and if h is big enough than the corresponding billiard 
map has non- vanishing Lyapunov exponents almost everywhere; moreover h —* oo as a —* y/2. In a 
previous paper we found a bound for h assuring the K-property for these billiards, for values of a 
■ very close to 1. 

In this work we study the stability of a particular family of periodic orbits obtaining a new bound 
for the chaotic zone for any value of a < y/2. 

o 
> 

X 

The elliptical stadium is a plane region bounded by a curve T, constructed by joining two half-ellipses, 
with major axes a > 1 and minor axes b — 1, by two straight segments of equal length 2h (see fig.0). 



1 Introduction 




Figure 1: The elliptical stadium. 
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The billiard on the elliptical stadium consists in the study of the free motion of a point particle inside 
the stadium, being reflected elastically at the impacts with T. Since the motion is free inside T, it 
is determined either by two consecutive points of reflection at T or by the point of reflection and the 
direction of motion immediately after each collision. 

Let s £ [0, L) be the arclength parameter for L and the direction of motion be given by the angle (3 
with the normal to the boundary at the impact point. The billiard defines a map T from the annulus 
A = [0,L) x (— 7r/2,7r/2) into itself. Let (so,A)) and (si,/?i) € A be such that T(s ,A)) = ( s iiPi) an d 
that T is C°° in some neighborhood of sq and si (notice that F is globally C 1 but not C 2 and piecewise 
C°°). Then, T is a C°°-diffeomorphism in some neighborhoods of (so,/?o) an d (si , It also preserves 
the measure g?/i = cos (5 dp ds (see, for instance, ||). 

{A, /x, T) defines a discrete dynamical system, whose orbits are given by 

O(s , fa) = {(s n , p n ) = r i (s , A), neZ}cA 

So the elliptical stadium billiard defines, almost everywhere, a two-parameter family of diffeomorphisms 
T a ^h whose dynamics depend on the values of a and h. For instance, when h = 0, T a j t is integrable for 
every a since we have the elliptical billiard. 

In p|, V. Donnay proved that the elliptical billiard stadium map T a ^ has non- vanishing Lyapunov 
exponents almost everywhere if 1 < a < \/2 and h is sufficiently large. He also proved that h must go to 
infinity as a approaches V2. Donnay addressed a challenge: "One could try to calculate bounds on these 
lengths." 

In 0], we proved that if 1 < a < \/ 4 — 2\/2, then h > 2a 2 V a 2 — 1 assures not only the positiveness of a 
Lyapunov exponent, but also ergodicity and the K-property. However, 2a 2 V a 2 — 1 does not seem to be 
an optimal lower bound for h. Numerical simulations exhibit chaotic phase spaces for values of h smaller 
than this bound. Moreover, v4 — 2y/2 ps 1.082 is far from a/2~ and close to 1, so in this case we are very 
close to the Bunimovich stadium which is chaotic for all h > 0. 

In this work, expanding and revising the work in we study a family of periodic orbits, with any pair 
period p > 4, whose behavior looks generic. This means that in phase space, when they are elliptic, they 
are surrounded by invariant curves which constitute elliptic islands of positive measure which disappear 
as the orbits change from elliptic to hyperbolic. 

We will prove the existence of a curve H = H(a), which diverges as a approaches V2 such that above it 
all these periodic orbits are hyperbolic. Even though H = H(a) may not be the optimal searched bound 
for chaoticity, it is at least a lower bound and seems not far from it. It also gives a very good answer to 
Donnay's challenge. 

This paper continues in the following way: in section 2 we describe the family of periodic orbits. In 
section 3, we study their hyperbolicity. Section 4 consists on the definition of the bound H(a). Section 5 
contains some concluding remarks. 



2 Pantographic Orbits 

On M 2 , we fix the origin on the center of the elliptical stadium and take the x-axes containing the major 
half-axis of the half-ellipses (see figure Q). 

Given a, h and a positive integer i, an (i, a, /i)-pantographic orbit, denoted by Pan{i, a, h), is a (4 + 2i)- 
periodic orbit, symmetric with respect to the coordinate axes, alternating the impacts up and down, with 
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exactly 4 impacts at the half-ellipses (2 at each one, joined by a vertical path) and 2i impacts at the 
straight parts (i at each one) and crossing the y-axis only twice (see figure ^) . 




i = 3 i = 4 i = 5 



Figure 2: Pantographic orbits of the elliptical stadium. 

The choice of those pantographic orbits was motivated by several remarks. The elliptical stadium billiard 
can be viewed as a perturbation of the elliptical billiard. The orbits of the elliptical billiard may be 
classified according to two different main features: those that have an elliptical caustic and those with 
a hyperbolic caustic (on the phase space of the elliptical billiard, these last orbits belong to invariant 
curves surrounding the elliptic period-2 point) . The invariant curves associated to the orbits with elliptical 
caustic are easily destroyed by the perturbation (in the same way it happens in the Bunimovich stadium) . 
This is not so easy for those with hyperbolic caustic. So orbits that have a chance to remain elliptic, 
after perturbation, must be close to those trajectories with hyperbolic caustic. 

On the other hand, it is known ]3j that the elliptic character of an orbit can be given by a relation 
between the curvature of the boundary at the impact points and the total length of the trajectory. Since 
bounces on the straight parts of the elliptical stadium only change the length of the trajectory, the elliptic 
behavior will depend fundamentally on the number of impacts with the elliptical parts. But, if a < no 
trajectory with hyperbolic caustic, on the elliptical billiard, can have more than two consecutive bounces 
on the same half-ellipse [|. Since while bouncing on the half-ellipse, any trajectory on the elliptical 
stadium billiard behaves exactly as a trajectory on the elliptical table, we must look for periodic orbits, 
close to orbits of the elliptical billiard with hyperbolic caustic and bouncing twice at the elliptical part 
of the stadium. 

Between all the periodic orbits with this behavior, the pantographic have the following properties: 

a) they exist for every even period p > 4 and so can be studied as the period goes to infinity. 

b) they can be explicitly localized and linearized. 

Strengthening our choice of this family, in the numerical simulations we have carried out, those panto- 
graphic orbits appear as the last ones having observable K AM- like-islands. 

Proposition 1 Pan(i, a, h) exists for every 1 < a < \[2 and h > 0. 
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Proof: Let P be the point of Pan(i, a, h) located in the right half-ellipse and on the first quadrant. Let 
A e [0, 7r/2] be such that P = (a cos A + h, sin A), and let j3 > be the angle of the trajectory by P, 
with the normal to the boundary. Using the obvious symmetries, it is easy to see that P is a point of 
Pan(i,a,h) if the straight line passing by P with slope tan(7r/2 — 2(3) cuts the y-axis at (0, — i). (see 
figure ||) 
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tan 2/3 



Figure 3: 

h + a cos A 



i + sin A 



And since tan (3 = cos A/ (a sin A) we have 



a tan A 



1 



a tan A 



2(n/l + tan 2 A + tan A) 
hvl + tan 2 A + a 



(1) 



It follows from the same arguments that a trajectory containing the vertical piece from P' = {a cos A + 
h, — sin X) to P — (a cos X + h, sin A) will cut the y-axis at (0, —y), where 



v = y(t) 



2+2 



=L=(Wl + t 2 + a) 



2aWl + 1 2 

with t — tan A. Then finding a solution of (|l|) is equivalent to find a t such that y(t) = i. 
It is not difficult to verify that lim t _>o y{t) — —oo, lim^oo y(t) = +oo and 
dy _ (1 + a 2 t 2 ) (a + h{l + t 2 ) VT+¥) 



(2) 



dt 



2at 2 (l + t 2 y- 



> for all fixed a > and h > 



Then, as y(t) is a continous strictly increasing function running from — oo to +oo as t runs from to oo, 
y(t) = i has a unique solution U(a, h) for each integer i, for every given a > 0, h > and so ([!]) has also 
a unique solution tj(a, h) = tan(Ai(a, h)). 
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On the other hand, as we know Q that for a < y/2 no trajectory crossing the x — axis between any 
two consecutive hits with the boundary can have three consecutive impacts on the same half-ellipse, we 
conclude that the next impact after P' and P on the vertical piece of the trajectory as described above 
must be on the straight part of the billiard. 

Then if s(Ai) is the arclength corresponding to the point P — (h + a cos Aj,sinAj) of the stadium, and 
< /3(Aj) = arctanl/(aij) < it/2 the orbit of (s(Aj), /3(Aj)) under the billiard map T is Pan(i,a,h). So 
is the orbit of (s(Ai),7r — /3(Aj)j. □ 

Remark: Using the same ideas, the existence of those Pan(i, a, h) can be proved for 1 < a < 2, but we 
are only interested on a < \/2. 

Lemma 2 Given i, a and h, let Xi = Aj(a, h) be the solution of Then A, goes to arctanl/a and 
/3(Aj) goes to ir/A as h goes to oo; Aj goes to it/2 as h goes to for all 1 < a < y/2. 

Proof: The right side of equation (|j) goes to zero as h — > oo so we must have in this limit, a tan A = 
l/(atanA) and it follows that tan A — > 1/a for all i. 

To study the behavior when /i — > it is sufficient to study for i = 0, because it follows from the proof of 
Proposition 1 that if i < j then ti < tj. For i = 0, equation ([!]) can be rewritten as: 

- (a 2 - cot 2 A) VT+ cot 2 A = [(2-a 2 ) + cot 2 A]cotA (3) 
a 

and when a < V2, h — > 0, [(2 - a 2 ) + cot 2 A] cot A — > and ij(0) -> 7r/2. Thus, U -> 7r/2 as /i -> 0. □ 

Remark: Let us call pantographic-like orbits on the elliptical billiard the periodic trajectories that 
have vertical segments both at left and right extremes. It would be amazing to compare the results in 
Proposition 1 with the existence of those pantographic-like orbits. As can be seen in the 2n-periodic 
pantographic-like orbit exists if a > a n where a n satisfies tan7r/n = 2-\/a 2 — l/(a 2 — 2). For instance, 
there is no 4-periodic pantographic-like orbit if a < y2, or 6-periodic if a < 2. 



3 Hyperbolicity of the Pantographic Orbits 



Proposition 3 For each i, let on = 




1. Fori > if a% < a < v2, there exists a unique hi(a) such that if h < hi(a), Pan(i,a,h) is elliptic 
and if h > hi(a), Pan(i,a,h) is hyperbolic. 

2. For i > 1, if I < a < cti , then Pan{i, a, h) is hyperbolic for all h > 0. 

Proof: For fixed a < \/2, i and h, let Aj(a, h) be the solution of (Q), (3 be the angle, with the normal, of the 
outgoing trajectory at P — (acosAi + /i,sinAj) and s the corresponding arclength. So T 4+2 ' (s, /?) = (s,{3) 
and to study the stability of this orbit we must analyze the eigenvalues of 2?T 4+2l |( s / 3). 

Let (s n ,f3 n ) and (s„+i, (3 n +i) be two consecutive impacts of a trajectory with the two different half-ellipses 
(with k > impacts with the straight parts between them), or two consecutive impacts of a trajectory 
with the same half-ellipse (with k = 0), then (see, for instance, |3)) 

pyfe+li _ (~l) fc f ln,n+lKn — COS /?„ l n ,n+l \ 

\(sn,0n) — CQS p n+1 y K n K n+ x{l n>n+ i - COS/3 n - COS^n+l) l n ,n+lK n +i - COS/3 n+ l / 
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where K stands for the curvature, and l n , n +i is the total length of the trajectory between the two impacts 
with the half-ellipses. 

Then, using elementary geometry and the symmetries of the trajectory, we can write 2?T 4+2l |( Sj/ 3) = 
(M 1 M 2 ) 2 with 



71/ 



3 cos (3\K 2 - 2 cos 0) l 3 K - cos /3 



1 ( lj K — cos (3 I j 



and where l\ = 2sinAi, l 2 = 2^(h + a cos A^) 2 + (i + sin A^) 2 and K = a/(a 2 sin 2 A^ + cos 2 \) 3 ^ 2 . 
Now if we define 

\cos/3 J \cos/3 ) 

we have that Pan(i, a, h) is elliptic if < Ai(a,h) < 1, parabolic if Aj(a, h) = or 1 and hyperbolic if 
Ai(a,h) < or Aj(a, /i) > 1, which means that the eigenvalues of PT 4+24 |( s /5 ) are respectively purely 
imaginary and unitary, equal to 1, real and one bigger than 1 and the other smaller than 1 (remember 
that the system is conservative). 

To study the function Ai(a, h), we need the following lemma: 

Lemma 4 The function Aj(a, h) has the following properties for h > and 1 < a < y/2: 



1. Ai(a,h) > 

3. lim/j^+oo Ai(a, h) = +00 



I lim^o Ai(a, h) = ^ (0) = - - 1) > 



Proof of the lemma: If a < y2 the half-osculating circles of the ellipse are entirely contained inside 
the ellipse 0], and so Zi > 2^£. Since h > h, A, is the product of two positive factors and property 1 
follows. 

To prove property 2, we remark first that, from formula ^ we can derive implicitly dt/dh and show that 
it is negative. 

We have 

hK 2 n l + t 2 



= 2- 



cos (3 a 2 sin 2 /3 + cos 2 (3 1 + a 2 t 2 
and this is a decreasing function of t if a > 1. So the first factor of Aj is a decreasing function of i. 

For 0<A<7r/2(i>0)ii'isa decreasing function of A, and so a decreasing function of t. As tan/3 = 1/at, 
cos /3 is an increasing function of t. This implies that Kj cos /? decreases with t. Now 

<9Z 2 —a sin A , 

-f = , tan 2(3 - tan/3 

aA 1 + sin A 

is negative, as [3 < n /A. 

Thus, Aj is a product of two decreasing functions of t, and so, is an increasing function of h which is 
property 2. 
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Property 3 is obvious since li — > oo as h — > oo and all the other quantities are bounded. 

When h — > 0, Ai — ► 7r/2, — > 2, cos/3 — > 1, K — > 1/a 2 and Z2 — * 2(1 + i), which implies property 4. □ 

Now we finish the proof of Proposition ^. Given i, we know from the lemma above that for each 
1 < a < a/2, Ai(a,h) is an increasing function of h running from Li (a) to 00. The function Lj(a) 
decreases with a, = l + 2i>lifi>0 and Lj(-\/2) = for every i. Li = 1 has the solution 

So, if cti < a < t/2) there exists a unique /ii(o) such that A,(a, /i;(a)) = 1 and /i < hi(a) implies 
Aj(a, hi(a)) < 1, h > hi(a) implies A,-(a, hi{a)) > 1. 

On the other hand, if a < a;, A, (a, h) > 1 for all /i and the result follows. □ 

For each i > fixed, and for all o>i < a < a/2, formula (|l|) and A,(a, h) = 1 constitute a system equivalent 
to 

a 2 sin 3 A + ^(a 2 -l)sin 2 A-sinA-^ =0 (4) 



hi ( a ) = si f A _ (2i sin A + 1 — (a 2 — 1) sin 2 A) (5) 

(a 2 + 1) sin A — 1 



where A = Aj(a). 



The values sinA = 1/a and h — y/a? — 1 satisfy the equations above for i = 0. So, /io(a) = Vr-T for 
1 < a < y/2. 

Now, let be y = sin Ai (a) . The problem of finding hi (a) is reduced to finding a root of the cubic polynomial 
in the interval (0, 1). Pi. a (y) — can be rewritten as 

*(» 2 -^) = i (1 - (a2 - 1)s2) 

The left hand side is a cubic polinomial with roots at 0, ±l/a; it is positive for y > 1/a and negative in 
(0, 1/a). The right hand side is a quadratic polynomial, with roots ±l/\/a 2 — 1 and which is negative for 
y > l/Va 2 — 1 > 1/a. This implies that Pi^ a has only one positive real root and that this root belongs to 
(1/a, 1/Va 2 — 1)- As the left hand side is for y = 1/a and the right hand side is positive, Pi )0 (l/a) < 0. 
However, Pi >a (l) > for a > on and Pj i0 (l) < for a > a,. This implies that for each i, Pi, a {y) has one 
and only one real root in (0, 1) for on < a < 

This root can be found by standard techniques: 

y t {a) = 2 ./A cosO - l{a l~ l) (6) 
6 a z 

where 

12a 2 + (a 2 -l)^ 2 



.4 
B 



(6a 2 ) 2 

2 (18 a 2 (l + 2 a 2 ) i + (a 2 - l) 3 i 3 ) 



(6a 2 ) 3 

cos 36» = -B/{2A 3/2 ) 
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Moreover, a more careful investigation shows that B/(2A Z / 2 ) < 1. When B/(2A 3 ^ 2 ) < — 1, there is only 
one real root and one should make use of the definition of cosine for imaginary arguments so "cos" is 
changed into "csh" . For — 1 < B/(2A 3 / 2 ) < 1, there are 3 real roots and we choose < < n/3, in order 
to have yi(a) > 0. 



One can also write (||) as 



A (-1 + a* 

" (o) = CF + C " "erf 



where C = (V B 2 — 4A 3 — B)/2, and when it is imaginary, the choice of the logarithm branch is the same 
as the choice of 9 above. 

If we introduce this value of sinAi(a) in (||), we obtain explicit formulae for the hi{a). These functions 
are plotted in figure ^. 




h 

A 



i = 1000 




Figure 4: Graphs of hi 



4 A lower bound for the chaotic zone 



For i = 0, ho(a) = Va 2 — 1 is a strictly increasing function in (1, \/2). Figure [| shows that this is also 
true for i = 1, 2. For i > 3 we have the following: 



Proposition 5 For each fixed i > 3, ^ > 0. 
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Proof: For i > 0, we have 



2ayy/l-y 2 /, -l + a 2 y 2 
-1 + (1 + a 2 ) y 2 1 



where y = y^(a) is the only root of Pj, a (j/) in (1/a, 1). Then, we have that 

dy -ay 2 (i + 2y) 



da i y (a 2 — 1) + 3 a 2 y 2 — 1 

Now 



< 0. 



can be written as 



As 



dhi dhi dhi dy 
da da dy da 



dhi dhi ( ^ a dy\ | f dhi a dhi\ dy 



da da \ y(l — y 2 ) da J \ dy y(l — y 2 ) da J da 



2yVT~^(l-y 2 -2a 2 y 2 + 3a 2 y 4 + a 4 y 4 - i 2 (I - y 2 + a 2 y 2 )) 
da ' " ^ z(-l + y 2 +a 2 y 2 ) 2 

dh t 2 a (1 - y 2 - 2 a 2 y 2 + 3 a 2 y 4 + a 4 y 4 + 2 a 2 y 4 - 2 a 2 y 6 - 2 a 4 y 6 - i 2 (1 - y 2 + a 2 y 2 )) 

~dy~ ' i VT^ 2 '(-1 + y 2 + a 2 y 2 ) 2 

we obtain that 

dhi a dhi — 4a 3 y 4 

dy~ y(l - y 2 ) ~da~ ~ \ JY~y^ (-1 + y 2 + a 2 y 2 ) 

and 

_a dy = 1 + y 2 (3 a 2 y 2 - l - a 2 ) + j y (1 - y 2 + a 2 y 2 ) 

y(l-y 2 )aa - (1 - y 2 ) (i y (a 2 - 1) + 3 a 2 y 2 - 1) 

Now we will show that dhi j da < for i > 3, so dhi /da > for i > 3. 



where 



^ = _*vVE}F( D - i *E) 

da i(-l + y 2 +a 2 y 2 ) 2 V 7 



D = (l-y 2 ) + a 2 y 2 (3y 2 -l) + a 2 y 2 (a 2 y 2 -l) >0 
E = l + (a 2 -l)y 2 >0 



Now 

S > 1 + ^ > 1 + ^ = ^ > 
3 6 6 

as V2 > a > ai = 2/V3, for i > 1, and 1 > y > 1/a > l/y/2. 
On the other hand 

< D < l - + 2a 2 y 2 + a 2 y 2 = - + 3a 2 y 2 < - + 6 = — 
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So if 



/39 D 



dhi/da < and the result follows. □ 
We were also able to found the asymptotical behavior of the hi at a = \ 2. 



Proposition 6 If a= \/2 then lim 



Proof: For a = \[2 we have 

% i % 

p i,V2(y) = y 3 + - 2 y ~ i 

so its only positive real root goes to 1 as i goes to oo. It follows that 

lim sin Ai(V2) = 1 . 



Since P i ^(sin A) = 0, 



sin A f sin 2 A | = - cos 2 A . 

2/4 



The left hand side goes to 1/2 as i goes to oo implying 



lim - cos 2 Ai(V2) = 1 . 

As 

, 9 / /s 2 cos 2 A . . . , 9 , \ 2 

h?( V2) = o (2* sin A + cos 2 A) 

1 ' 3sin 2 A-l V ' 

using the limits above, gives the desired result. □ 

Proposition ^ tells us that there exists an N such that if i > N, hf(y/2)/4 sa i. Plotting hf(^/2)/4 x i, 
we have observed the same linear behavior also for small values of i. 

For each 1 < a < \J~2 there exists a j such that a.j-% < a < atj. So, Pan(i,a,h) is hyperbolic for i > j 
and is elliptic if < i < j and h < hi{a). There exists, then, only a finite number of h^s defined for 
this value of a and we can define the announced lower bound by H(a) = maxj<j {hi(a)}. Proposition [j] 
implies that H(a) — > oo as a — * \/2. H{a) can be seen on figure 0. 



5 Final remarks 

Clearly above the curve H (a) all pantographic orbits are hyperbolic and bellow it some have eigenvalues 
in the unit complex circle. Generically, KAM theory establishes the existence of positive measure elliptic 
islands surrounding those orbits and our numerical simulations corroborate this result. In this paper we 
have not proved the existence of such islands. 

On the other hand, above this bound, although we can prove that all the orbits of this family are 
hyperbolic, we can not assure the non existence of other elliptic periodic orbits that could be surrounded 
by positive measure sets of invariant curves and having vanishing Lyapunov exponents. As a matter of 
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fact, we have observed other periodic orbits of T a ,h but the elliptic islands around them seem to disappear 
for values of h < H(a). 

A numerical case study which seems generic for this problem is presented in figures §, | and 0. There, 
we show the phase space associated for a — 1.24 fixed and different values of h. For this value of a we 
have that a 2 = Vl-5 < 1.24 < y/1.6 = a% (see fig. [|) so Pan(i, 1.24, h) are hyperbolic for any h > and 
i > 3. The pantographic orbits of period 4, 6 and 8 (Pon(0, 1.24, h), Pan(l, 1.24, h) and Pan(2, 1.24, h)) 
are the only relevant pantographic orbits as they are elliptic for small h. In fact we have ho ~ 0.7332, 
hi « 1.0236, h 2 w 0.6770, h 2 <h <h x = #(1.24). 

Note that for very small h (in this example, 0.1), a very rich structure of elliptic islands (the white holes 
on fig. can be observed; some of these islands correspond to other periodic orbits. 



Figure 5: a = 1.24, h = 0.1. 150,000 iterations of a single initial condition. 

As h is increased (fig. ^), they gradually disappear and only the pantographic islands seem to remain 
at h = 0.45. Then, at h = 0.73 > h 2 , the eight islands around Pan(2, 1.24, h) have disappeared and at 
ho < h = 0.75 < hi only the six islands around Pan{\, 1.24, h) can be seen, since Pan(0, 1.24, h) has 
also become hyperbolic. In each case, the region outside the elliptic islands seems to be a single ergodic 
component as it is filled up by a single orbit. 

For h = 1.05 > hi = H{1.2A) (fig. ^) the system seems to be ergodic. 

As we can see on figure ^, the order of extinction of the pantographic elliptic islands depends on the value 
of a. As shown above, for a = 1.24, the order is i = 2, then i — and then i = 1. But, for values of a a 
little bigger than 1.24, it can be 0,2,1 or even 0,1 and 2. 

However one can learn from figure ||, that as a approaches \[2 the last islands to disappear correspond 
to orbits with long period. This is expressed in proposition ^| So in this system, and this fact seems 
also to occur in other problems, the obstruction to ergodicity seems to be existence of elliptic islands 
around orbits with arbitrarily large period and so the existence of an arbitrarily large number of very 
small islands. Those islands, even though summing up into a positive measure region, can be invisible in 
a numerical simulation. 

We would like to point out that the pantographic orbits also seem to be the most important ones from 
the point of view of their focusing properties. In H, the existence of a positive Lyapunov exponent 
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Figure 6: a = 1.24 and different /i's. On the left side, 150,000 iterations of a single initial condition and. 
On the right, the iteration of a few initial conditions close to the elliptic pantographic orbits. 
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Figure 7: a = 1.24, h = 1.05. 150,000 iterations of a single initial condition. 



for the elliptical billiard was obtained through the study of the behavior of the caustic pencil, or the 
tangent vector to the invariant curves of the elliptical billiard. In the pantographic orbits, after hitting 
a half-ellipse twice, in the vertical portion of the trajectory, it crosses to the other half-ellipse. At this 
moment the focusing distance of the caustic pencil can be very large. This lack of focalization should 
be compensated by a larger traveling distance, so a bigger h, in order to have a splitting of neighboring 
trajectories. However, as a approaches \f2 this focusing distance may tend to oo. The loss of ellipticity 
by all the pantographic orbits caused by increasing h indicates that the behavior of the caustic pencil may 
be controlled at this point and one should be able to prove the existence of positive Lyapunov exponents. 
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